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Page U; In the second line of equation ( 2 ), the quantity Ayy^ which 
appears in the numerator of the fraction under the radical in the first 
term within the brackets should be AZy^. 

Page 6 : In equation ( 8 ), the minus sign inside the parentheses should be 

a plus sign. 


Page 6 : In equation { 9 ) > the minus sign inside the parentheses should be 

a plus sign. 




4 


t 

f 


r 

h‘ 



y 





( 


L 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2553 

PITCHING-MOMENT DERIVATIVES Cm AND C^^ AT SUPERSONIC 
SPEEDS FOR A SLENDER-DELTA-WING AND SLENDER-BODY 
COMBINATION AND APPROXIMATE SOLUTIONS FOR 
BROAD-DELTA-WING AND SLENDER- 
BODY COMBINATIONS 
By Arth\ir Henderson, Jr. 


SUMMARY 

The pitching-moment derivatives Cm^ arid C^^ at supersonic speeds 

are developed for a slender -delta-wing and slender-hody combination 
having no afterbody. By drawing an analogy between the aerodynamics of 
the wing-body section of the combination and the aerodynamics of a delta 
wing alone, the results for the slender-delta-wing and slender-body com- 
bination are modified to the extent that approximate solutions for C^ 

and Cm^ for broad-delta-wing eind slender-body combinations can be 

obtained. 


INTRODUCTION 


Various methods, based on linear theory, for obtaining solutions 
for the flow about wing-body combinations have been developed for the 
determination of the lift and moment due to angle of attack. Refer- 
ences 1 to 7 comprise a fairly comprehensive list of most of the signif- 
icant of these methods, which include both approximate and exact solu- 
tions. All the exact solutions to the linearized differential equation 
of steady supersonic flow, however, employ iteration processes, infinite 
series, or both, and their practical application results in approximate 
solutions although the error is often negligible, depending upon the 
particular problem, rate of convergence, number of iterations, and so 
forth. Spreiter (reference 7) has presented solutions in closed form 
to the two-dimensional Laplace equation of potential flow for the lift 
and moment of wing -body combinations. These solutions apply to the super- 
sonic range for the limiting case of a slender wing-body configuration. 
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For the stability derivatives of wing-body combinations, there are 
a few papers on the damping-in-roll characteristics (see, for example, 
references 8 and 9) but none for the damping in pitch. 

The purpose of the present paper is to extend the method used by 
Spreiter in reference 7 to the calculation of the pitching-moment deriva- 
tives due to constant rate of pitch ^^mq\ to constant accel- 

erated motion in the vertical direction (Cny^] for a slender-delta-wing 
and slender -body combination. In addition, an approximate solution to 
these derivatives is developed for a broad-delta-wing and slender-body 
combination in supersonic flow by introducing certain modifying factors 
into the slender-delta-wing and slender-body resiilts. 

Certain conditions are placed upon the configvLration. The body 
ahead of the wing is slender, has a circular cross section, and is 
pointed at the nose, and the slope of the body meridian section is 
continuous. For the wing -body section, the wing semiapex angle is small) 
along the wing-body juncture, the body radius is a maximum and is con- 
stant; and finally, the configuration has no afterbody (see fig. l) . 


SYMBOLS 


0 , 0 ' 

Z 

R 


a 

s 

w 

r,0 

q 


potential functions 
stream functions 
complex variable (y + iz) 

body radius (R = R(x) on body ahead of wing and R = a along 
wing-body section) 

body radius along wing-body section 

y-coordinate of wing leading edge 

velocity in positive z-dlrection 

polar coordinates 

constant angular velocity of pitch 
constant time rate of change of angle of attack 
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p 

p 

t 


x,y,z 

V 


n 

M 

A 


C 


m 


C 


niq 


C 


“Id 


c 

c 


c ' 
I 


perturbation pressure (difference in pressxxre between body- 
surface and free stream) 

density of fluid 

time 

Cartesian coordinates 

free -stream velocity 

point of rotation measured from nose 

inward- drawn unit normal -vector 

pitching moment 

area of basic wing (including portion enclosed by body) 
pitchlng-momerrt coefficient 




nondimensionaJL stability derivative due to constant rate of 



nondimensional stability derivative due to constant accelerated 
motion in vertical direction 

root chord of basic wing 
mean aerodynamic chord of basic wing 

root chord of exposed wing 
total length of wing-body configuration 
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semiapex angle of basic wing 
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max 


h,m 

M 


value of s at x = Xq or at Cq 

maximum value of s (value of s at x = 2) 

point of rotation measured from apex of basic wing) positive 
in positive x-direction 


interference factors 
Mach niomber 


p = _ 1 

K = tan e 

Q constant of integration 


k = 


max 


E'(PK) complete elliptic integral of second kind 


F'(pK) 


X-i = 


(I 


sin^Q d0 


complete elliptic integral of first kind 

>„/2 d 0 


f rn/ 2 \ 

I |A - (l - 3%^)sin2ej 


1 - 


(l - 2P%2^E'(PK) + P%^'(pK) 


^2 E'(PK) 


3 p2 2 ^2 1 
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Subscripts: 

W wing 

B body 

q due to q 

d due to d 


AMLYSIS 


The linearized differential equation of steady supersonic flow is 


2 ^^0 ^^0 1 ^0 1 




3r 2 




= 0 


At present an exact solution to this equation does not exist in closed 

form for wing -body combinations. However^ if the term ^ becomes 

very small with respect to the other terms of this equation^ it may be 
neglected. Solutions to the Laplace equation which results from dropping 

2 

the term P — ^ have been found in closed form for the lift and moment 
due to angle of attack (reference 'J) . It has been found that the condl- 

tion necessary for P^ ^ to be negligible for the angle -of -attack case 
is that the configuration be slender and that P be not excessive. For 


a delta-wing and body combination^ the term slender implies that 

d% 


6R 

dx^ 


dx- 


— ^ and K are very small. 


In the present paper ^ which treats the steady-pitching and the 
time-dependent^ constant -acceleration cases of delta-wing and body 
combinations^ a velocity potential satisfying the two-dimensional 
Laplace equation is used. In the appendix it is shown that the conditions 
to be satisfied for the Laplace solution to be applicable to the super- 

2 

sonic range are that ^ very small. 

dx^ 
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After a velocity potential which, satisfies the Laplace equation is 
found, the next step in the analysis is the determination of the pressure 
distributions over the slender -de It a-wing and slender-body combination 
restating from the two types of motion which give rise to and 

namely, constant rate of pitch and constant accelerated motion in the 
positive z-direction, respectively. When the pressure distribution is 
known, the moment may be calculated about any axis of the conflgtu’ation, 
and, from their respective definitions, Cjjj and Cjjj^ may then be deter- 


mined. The configuration to be considered and the coordinate system 
employed are shown in figure 1. 


I 


Velocity Potential 


Spreiter (reference 7) shows that the complex potential for a 
\miform stream of velocity w at infinity flowing vertically downward 
over a stationary two-dimensional circular cylinder symmetrically located 
on a horizontal flat plate is 


0 ’ + 




( 1 ) 


where 


Z = y + Iz 

R radius of cylinder 

s semispan of plate measured from center of cylinder 


For a slender configuration describing a slow, steady pitching 
motion, the cross-flow velocity distribution is, to the first order, 
proportional to x. Inasmuch as potential flow is assumed, this velocity 
distribution must be looked upon as being generated by the motion of the 
configuration ,ln fluid which is at rest, because, if the distribution 
were due to the motion of the fluid about a stationary body, the flow 
must be rotational and the assiimption of potential flow is then violated. 


The complex potential of the aforementioned configuration moving 
upward through still air with the vertical velocity w then is 



0 + i\|f = iw- 


( 2 ) 
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Transforming to polar coordinates (Z = r(cos 6 + ± sin ©)) and solving 
for the velocity potential gives 


Equation (3) is the general expression for the velocity potential. 

Whether 0 pertains to the constant-pitching or the constant -acceleration 
case depends upon the value of w. For a wing-body configuration pitching 
about a point from the nose, the vertical velocity w varies along 

the length of the configuration according to w = q (x - Xq). For con- 
stant acceleration in the positive z-dlrectlon, the velocity varies with 
time according to w = dVt. 



1/2 



(3) 


Pressure Distribution 


The equation for the pressure distribution is 



The term 



does not contribute to either the lift or moment since 
symmetric and on the wings, although ^ is antisym- 


on the body it is symmetric and on 



is symmetric; therefore, for the configuration considered. 



( 5 ) 


For the case of pitching with constant angular velocity, 
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and for constant acceleration, evaluated at time t = 0, 


PI = P 


"a 




(7) 


In order to determine the loading over the wing-body combination as 
given by equations (6) and (7)^ the press\rre distributions in two regions 
must be considered for each expression. They are: 

(a) Pq^a on the body where r = R 

P. • on the wing where 6=0 and a "S r = s 
q,(x 

For the pitching wing -body combination with w = q(x - Xq) and the 
preceding conditions, equation (6) gives for the pressure over the body 
and wing, respectively. 



= -pVq 


-Rs sin 6 + ' 

/s^ + R^ - 2r2s2cos 26 


s 


+ 


1 

X 

( r 2 - s2cos 26 ) 

^ dR , (x - xq)I 

[ 8^ - R^: 

1 ds 



sl4^ + R^ - 2 R^s2cos 26 s^l/a^ + R^ - 2R2s2cos 26 


(8a) 



2r3(x - Xp)(r^ - 8 ^) dR 
rs/(r2 - s2)(r^ - r2s2) 


(8b) 


Similarly, for the case of constant acceleration, in which w = Vat is 
used, equation (7) yields 


/ \ 1/4^ + R^ - 2R^s2cos 26 - Rs sin 6 

(p^)b = -PVa 


s 


(9a) 


,2Z 
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(w)„ - -pvi 


|/7~- 

rs 


(9b) 


Cj„ and for Slender -Delta-Wing and Slender-Body Combinations 

The moment on the wing-body configuration measiired about a point 
from the nose is 


M = / (x - X )n . p dA 


( 10 ) 


where n is an inward-drawn unit vector, normal to the surface, and A 
represents the surface area of the configuration. Now 


M 


^m 1 p - 

ipV^Ac 


anci 



Therefore Cmg and Cn,^ are respectively. 


C = ^ -ii-J 
“q ^IpVAc^ 


■'Z-c' P2jr 


JO Jo 


(x - Xo)(pq)BH sin 0 d0 dx + 


/Z-c ' Ja 


'Z PIT/ 2 

(x - Xo)(p ) 

-c'Jo ^ 

(x - x^)(pq)„ dr dx 


R sin 0 d0 dx + 


( 11 ) 
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Xo)(P(^)bR sin e d0 dx + 



(x - Xo)(p^)gR 


sin © dO dx + 



( 12 ) 


where the first integral in each expression is the contribution of the 
body ahead of the wing and the last two in each expression are the 
contributions of the wing-body section. 

The conditions to be imposed in evaluating these integrals are: 

(a) On the body ahead of the wing, s = R 

(b) On the body at the wing -body section, R = a and — = 0 

dx 

(c) On the wing, — = Constant = tan e 

dx 

Integration of the terms for the wing-body section of the configuration 

may be simplified by making the substitutions for x and Xq which are 

suggested by condition (c). Since — = tan e, 

dx 

s = X tan e + Q 

where Q is a constant. Therefore, 

X = ^ ~ 

tan 6 

X = ~ Q 

o tan e 


tan € 
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and the limits of Integration are now from s = a to s = s^ ^y . From 

the geometry of the configuration (see fig. l)^ %ax = ^ ^ and 

Sq “ ^o ^ ) where Cq Is the location of the point of rotation 
measured from the apex of the basic wing, positive In the positive 
x-dlrectlon. 


Performing the operations Indicated In equations (ll) and (12) and 
substituting limits, = c tan e, and Sq = c^ tan e resvilts in 




(x - Xo)R^dx 



(x - Xq)% — dx 
^ dx 


6n tan € ( ^ e 



+ kn tan e 




(13) 



(x - Xo)R^dx - 2it tan 




(1^) 


where 


e = 1 - 


3 ^ 


i^3 3 


loge 


f = 1 




g = 1 - 2k 


2 + k^ 


h = 1 - kk^ + (3 + ^ loge ijk^ 


m 


= 1 - 6k^ + 8k3 - 3 k^ 


The variation of these interference factors with k is shown in figure 2. 
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Equations ( 13 ) and (l4) are the expressions for and Cju^ for 

a slender-delta-wing and slender-body combination corresponding to the 
conditions stipiilated. When these terms are added to obtain the damplng- 
in-pltch parameter 0^^ + Cnj^, integration by parts allows the resulting 

expression to be written as 






6jt tan ^(0 ® 


4jt tan e gj - 2rt tan h - mj 


(15) 


Again, from the geometry of the configuration, when k / 1, 


(Z - c') - Xq = -(Cq - 


° tan ej 

This relation allows equation (I 5 ) to be written as 


Cjn + Cm, = -4nk'^ tan e 




I k2 - 3k 


c \ c , 


- 6jt tan ^ I ^ ^ 




4jt tan € -^(f - g I - 2it tan ^(0 H 


( 16 ) 


When k = 1, the wing span goes to zero, and for a slender body of 
revolution 


"mn = 


^ J (k - k,)H=ax - ^ jT' (X - x.,% g .X U7) 


"“i = ■ Ac 2 jo 


(x - XQ)R^dx 


( 18 ) 
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and 



+ C 




kna^, 
2 '' 


Ac 



(19) 


where A and c represent some characteristic area and lengthy re spec- 
tlvely, of the body. Equation (19) agrees with Miles' result (refer- 
ence lO) if A = Jta^ and c = Z. 

When R = a = 0, the body radius goes to zero, and for a slender 
delta wing 

Cmq = tan + krt tan e ^^1 - (20) 

Cn^ = -2n tan € (21) 

which are the expressions for Cmg and for the slender delta wing 

fmind by Ribner (reference ll) . 

From these equations for ^mg the terms for the wing- 

body section of a slender-delta-wing and slender-body combination are 
seen to be in the same form as Cjug and for the basic wing alone . 

Each term of the equations for the basic wing alone is modified by a 
factor which is a function of the ratio of the body diameter to the 
maximum wing span. This modification is due to the interference effects 
which result from placing a slender body on a slender delta wing. 


C™ and C . for Broad-Delta-Wing and Slender-Body Combinations 
“q ^ 

From practical considerations, solutions for *^mg 
broad-delta-wing and slender-body combinations in supersonic flow are 
desired. A method of obtaining an approximate solution to this 
problem from the preceding development is s\;ggested by the similarities 
between the expressions for the slender delta wing alone and for the 
slender-delta-wing and slender-body section of the configuration. An 
intuitive approach would be to assinne that a delta-wing and slender- 
body section, in going from a slender-delta-wing and slender-body section 
to a broad-delta-wlng and slender-body section, follows the same laws 
that a delta wing alone follows in making the same transition ( see the 
next section for a discussion of the validity of this assumption) . 
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Investigations by Brovm and Adams (reference 12) and by Ribner and 
Malvestuto (reference 13) made after the publication of Ribner 's paper 
on the stability derivatives of slender delta wings (reference 11) show 
that the stability derivatives of broad delta wings in compressible 
supersonic flow such that |3 tan e < 1 are the same as the results for 
the slender delta wing multiplied by certain elliptic integrals which 
are functions of the wing semiapex angle and the Mach number of the flow. 
Applying these laws to the wing -body section gives 



where ^2) ^3 the appropriate elliptic Integrals ( see 

fig. 3). The damping -in-pitch parameter is 



+ Cjj^ = -4nk tan 





tan 




tan 6 



+ tan 




(24) 


In order to determine approximate expressions for Cj„ and 

for the configuration when the wing leading edges are supersonic 
(P tan e > 1), the analogy drawn previously between the laws followed 
by a broadening delta wing alone and a broadening- delta- wing and slender- 
body section is continued into the region where P tan e > 1. 

As a delta wing alone continues to broaden to the extent that 

P tan f > 1, the equations for Cm (see reference 12) and Cm. 

q ^ 


are 
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was obtained by use of equation (I 5 ) in reference 13 and agrees 

with Miles' result (reference l4) Therefore the derivatives for a 

broad-delta-wing and slender-body combination in supersonic flow, such 
that 3 tan e > 1, may be approximated by 




XQ)R^dx 



(x - Xo)2r H dx - 




(27) 


and 

^ (x - Xo)R^dx + h - m^ (28) 

provided the body ahead of the wing-body section remains slender with 
respect to the Mach cone emanating from its nose. 


and 


Because of the nature of the factor 
Cm^ for P tan ^ 1 ^ a general curve^ 


and the values of 
such as Crn + Cr«. 



plotted against P tan cannot be drawn, 
have therefore been chosen and curves of Cj^ 


Certain basic delta wings 

+ Cjjj* plotted against 
q ^ 


M 

V 


have been drawn for different values of k. These curves are presented 
in figure 4. 
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DISCUSSION AND CONCLUDING EEMAEKS 


By an extension of the method used by Spreiter in reference 7 the 
pitching -moment derivatives C^j^ and Cjn^ for supersonic speeds have 

been developed for a slender-delta-wing and slender-body combination 
having no afterbody. By drawing an analogy between the aerodynamics of 
the wing -body section of the configuration and the aerodynamics of a 
delta wing alone^ the results for the slender-delta-wing and slender- 
body combination were modified to the extent that approximate solutions 
for Cjjj and Cjn^ for broad-delta-wing and slender-body combinations 

were also obtained. 


In order to check the validity of the reasoning used *in arriving at 
the assumption by which the approximate solutions were obtained, the same 
reasoning was applied to Spreiter *s results for the lift-curve slope Cj^ 

of a wing-body combination for which an exact solution to the linearized 
supersonic -flow equation also exists (reference 6). 


In reference 6, Browne, Friedman, and Hodes have presented an exact 
solution to the linearized equation of steady supersonic flow for a 
delta-wing and slender-conical -body combination for which the apexes are 
coincident. Spreiter (reference 7) tias presented a solution to the two- 
dimensional Laplace equation for the same configuration. In order to 
obtain some indication as to the reliability of the assumption made, the 
same reasoning was applied to Spreiter results for of the delta- 


wing and conical-body configuration as was applied to the C^, and C 




results of this paper, and the modification of Spreiter’ s results were 
then compared with the results of reference 6. The results of this 
comparison are shown in figure 5 wherein is plotted against 


P tan € for different values of k. For k = O. 7 O the curve from 
reference 6 is incorrect for high values of P tan € because an insuf- 
ficient number of terms of the series resiilts were taken. 


From the results of this comparison it appears that values of 
and for broad-delta-wing and slender-body combinations will give 

fairly good approximations up to at least k = O. 5 O. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., August 21, 1951 
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APPENDIX 


CONDITIONS FOR LAPLACE SOLUTION TO APPLY TO SUPERSONIC RANGE 


In the limit, as P 


2 ^^0 


2M 

V 




V2^t2 


approach zero, 


a solution to the two-dimensional Laplace equation 


^ 1 1 


is a solution to the linearized equation of supersonic flow 

2M^ 

V Sx 


Therefore these two equations are compatible if the above limiting con- 
ditions are satisfied. 

In equation (3) a solution to the two-dimensional Laplace equation 
is given as 

0 = 0(v,R, s,r,0) (29) 

where R = R(x) and s = s(x) = Kx + Q. If, for the present, the assump- 
tion is made that w = w(x,t), from- equation (29) 

^^0 _ ^ ^ 2 / ^^0 ^ ^ ^w ds I 

^^2 ^^2\Bxy ^^2\^xJ 3 g 2 \Bx/ \^w Sr Sx dx Sw Ss Sx dx 

S^0 dR ds \ ^ ^ ^ d^R ^ d^s 

Sr Ss dx dx I Sw ^^2 Sr ^^^^2 Ss ^^^2 

S^0 _ S^0 Sw Sw S^0 Sw dR S^0 Sw ds ^ ^ S^ 

Sx St ~ 3 y 2 Sx St Sw Sr St dx Sw Ss St dx Sw Sx St 


2 S ^0 S ^0 1 ^ _ 1_ ^^0 + 

Sx^ Sr^ ^ ^ r^ Sr^ 


0 ^ ^ 

St y2 ^^2 


(31) 
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and 

For constant rate of pitch, v = w(x) = q(x - Xq), and from 
equations (30) to ( 32 ) 


(32) 



at ^^2 




For constant accelerated motion in the vertical direction, 
w = w(t) = aVt, and from equations (30) to (32) 


j2 ^ . p2 


^ k2 + 2 ^ K = + 




2\dx 


dB . 30 d% 




dR os dx ^ . 2 
dx 




2 ^ 


a2( 32(Z 


dR 


V = 


32<Z 


oK 


y (34) 


M' 


2 ;.2, 


= m 2 ii a2 


3t^ 3 w2 


An examination of equations (33) (3^) shows that^ in order for 

the Laplace solution to be a solution to the linearized equation of 

2 

JTD * 

supersonic flow, — , ■= — , K, q, and a must approach zero. 

dx2 


Within the framework of the small-disturbance theory, however, such 
stringent conditions as these are not necessary for the Laplace solution 


dR d2R 



to apply to the supersonic range. Rather it is required that 
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K, aad a be of such an order of magnitude that P 

m 2 d^fl , ^ 2 , 

v2 at^ 

terms of the linearized equation of supersonic flow 


and ^ he negligibly small compared with the remaining 


I 
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Figure 1,- Geometry of configuration and coordinate system used. 



Figure 2.- Variation of interference factors vith k. 
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Figure 3«- ^3 and variation of and with p tan e. 
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(a) € = 15° • 

Figure U.- Variation of Cjjj + C_. with Mach number for various values 

(J *“(X 

of k for -^ = 0.85. 
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(b) € = 30°. 

Figure 4.- Continued. 
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(c) e = 45°. 
Figure 4.- Continued. 
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(d) 6 = 6o°. 

Figure 4.- Concluded. 
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Figure 5.- plotted against ^ tan e for various values of k for 

delta-wing and conical-body combination. Comparison of exact solution 
of reference 6 and modified solution of reference 7, 
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